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On the level of single atoms and photons, the coupling 
between atoms and the electromagnetic field is typically 
very weak. By employing a cavity to confine the field, the 
strength of this interaction can be increased many orders of 
magnitude to a point where it dominates over any dissipa- 
tive process. This strong-coupling regime of cavity quantum 
electrodynamics (QED)-^'^, has been reached for real atoms 
in optical cavities^, and for artificial atoms in circuit QED^ 
and quantum-dot systems^. A signature of strong coupling 
is the splitting of the cavity transmission peak into a pair of 
resolvable peaks when a single resonant atom is placed in- 
side the cavity - an effect known as vacuum Rabi splitting. 
The circuit QED architecture is ideally suited for going be- 
yond this linear response effect. Here, we show that increas- 
ing the drive power results in two unique nonlinear features 
in the transmitted heterodyne signal: the supersplitting of 
each vacuum Rabi peak into a doublet, and the appearance 
of additional peaks with the characteristic ^/n spacing of the 
Jaynes-Cummings ladder. These constitute direct evidence 
for the coupling between the quantized microwave field and 
the anharmonic spectrum of a superconducting qubit acting 
as an artificial atom. 

Circuit QED realizes the coupling between a super- 
conducting qubit and the microwave field inside an on- 
chip transmission line resonator^'^. The drastic reduction 
in mode volume for such a quasi- Id system^i^, the re- 
cent improvement in coherence times^^, and the absence 
of atomic motion render circuit QED an ideal system 
for studying the strong-coupling limit. In the custom- 
ary linear-response measurement of the transmitted mi- 
crowave radiation, the vacuum Rabi splitting manifests 
itself as an avoided crossing between qubit and resonator. 

It may be argued that the observation of the splitting 
is not yet a clear sign for quantum behaviou r^^i^^ , since 
avoided crossings are also ubiquitous in classical physics. 
However, quantum mechanics gives rise to a distinct non- 
linearity of these splittings when initializing the resonator 
in a higher photon Fock state: when the resonator mode 
is occupied by n photons, the splitting is enhanced by 
a factor \^n + 1 as compared to the vacuum Rabi situ- 
ation. This nonlinearity has been observed with single 
atoms in both microwave^^ and optical cavitie^i^. In 
circuit QED, this characteristic trait has been observed 
very recently both spectroscopically^^, as well as in time- 
domain measurement^!^. 

Here, we present a theoretical analysis and experimen- 
tal investigation of the vacuum Rabi resonance in a cir- 
cuit QED system, where we study the ^/n nonlinearity 
up to n = 5 by exploring the power dependence of the 
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Figure 1: Extended Jaynes— Cummings level diagram 
of the resonator— transmon system, a, Bare levels in the 
absence of coupling. The states are denoted |n, j) for photon 
number n and occupation of transmon level j. b, Spectrum of 
the system including the effects of transmon-resonator cou- 
pling. The effective two-level system relevant to describing 
the lower vacuum Rabi peak comprises the ground state and 
the antisymmetric combination of qubit and photon excita- 
tions. 



heterodyne transmission (see Methods). This type of de- 
tection is particularly simple, since it is a continuous mea- 
surement involving only a single driving frequency. We 
find that, surprisingly, each vacuum Rabi peak develops 
into a new doublet (supersplitting), which can be simply 
explained by the saturation of a two-level system com- 
prising both photon and qubit degrees of freedom. For 
even stronger driving, a series of additional peaks due 
to multiphoton transitions^^ up the Jaynes-Cummings 
ladder appear, but the complete shape of the spectrum 
closely matches the predictions from a master equation 
treatment including only relaxation, dephasing, and the 
known multilevel spectrum of the superconducting qubit. 
The quantitative agreement confirms the validity of the 
extended Jaynes-Cummings Hamiltonian for our system, 
rendering circuit QED an excellent tool for probing the 
fundamental interaction of matter and light. 

We consider a system consisting of a superconducting 
transmon qubit ^^^^^ whose multiple quantum levels are 
coupled to a single mode of a transmission-line resonator. 
Our arguments, however, will be general and should ap- 
ply to any other realization of cavity QED with suffi- 
ciently large coupling and atom coherence times. Taking 
into account the relatively small anharmonicity of the 



transmon qubit, the system is modelled in terms of an 
extended Jaynes-Cummings hamiltonian including a co- 
herent drive at frequency Cc;d/27r, 

H = hu;ra'^a^^hu;j\j){j\ (1) 

j 

+ Yl % (^ b' + 1) 01 + h.c.) + h^a^e-'^^' + h.c). 
j 

Here, the operator a {a^) annihilates (creates) a photon in 
the resonator mode with frequency uj^l2i\. The energies 
corresponding to the transmon eigenstates are hjjj {j = 
0,1,...), and they are determined by the Cooper pair 
box hamiltonian^ for given Josephson energy Ej and 
charging energy Ec- In the transmon regime considered 
here, Ej/Ec ^ 1, the levels form a weakly anharmonic 
ladder with anharmonicity a = 0012 — cjqi ^ —Ec {ooij = 
ujj — uoi is the transition frequency between levels i and 
j). The transmon-resonator coupling strengths are Qj ^ 
9oVJ + I5 where ^o/27r is the vacuum Rabi frequency. 
The resulting energy level diagram for vanishing drive, 
^ = 0, is schematically shown in Fig. [TJ 

The interaction of the cavity-transmon system with its 
environment enables control and measurement by means 
of microwave drive and detection^!^, but also leads to 
decoherence. The appropriate framework to capture all 
these effects is a description in terms of the reduced den- 
sity matrix p for the cavity and transmon. The Marko- 
vian master equation which governs the evolution of p is 
given by 



P = -^[H,p]^ KV[a]p + 7iP [Yl "^J \^) + 1| 
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where V[A] is the usual Lindblad damping superopera- 
tor defined by V[A]p = {[Ap,A^ + [A, pA^) /2. The 
three damping terms on the right-hand side model the 
loss of cavity photons at rate tz^ the intrinsic relaxation of 
transmon excitations at rate 71, and the pure dephasing 
of transmon state superpositions at rate 7^^ and relative 
strengths Ofj, Pj (see Methods). 

The quantity of interest for detection of the vacuum 
Rabi splitting is the transmission amplitude A{u;d) of 
a microwave signal at frequency uj^j^ii. Here, we em- 
ploy a heterodyne measurement of the field quadratures 
^ = (^out + ^Lt) and Q = (^ih\^^-ih^^^), where 6out refers 
to the output field at the drive frequency. Input-output 
theoryS^ enables us to express the steady-state transmis- 
sion amplitude as 



A = ^JP + Q2 = 2v^Ka)| = 2y^|tr(aps)|. (3) 

Our calculations thus only require the steady-state solu- 
tion ps of the master equation ([2|) . We explore the trans- 
mission intensity A^ as a function of the drive strength 
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Figure 2: Supersplitting of the vacuum Rabi resonance 
Avhen probing heterodyne transmission beyond linear 
response. The experimental data are obtained with a circuit 
QED system in the strong-coupling regime, where the vacuum 
Rabi splitting is observed to exceed 260 linewidths, see a. All 
plots show the squared heterodyne amplitude A? in arbitrary 
units, b, Measured transmission (colour scale) for the left vac- 
uum Rabi peak, as a function of drive frequency and power. 
The plot reveals the supersplitting of a single Lorentzian into 
a doublet of peaks, c— f, Cuts for constant power at the values 
indicated in b. In linear response, f, the vacuum Rabi peak 
is Lorentzian; as the power increases a central dip develops, 
e, leading to supersplitting of the peak, d, and eventually 
becoming asymmetric at the largest powers, c. The experi- 
mental data (red line) is in excellent agreement with theory 
(black line). The results from the 2-level approximation are 
shown for comparison (green dashed line). 



^, and compare our results to experimental data ob- 
tained for a transmon qubit in resonance with a 00^/27: = 
6.92 GHz coplanar waveguide (CPW) resonator with cou- 
pling strength Qq/ti = 347 MHz (see Methods for exper- 
imental details). In the linear-response regime, the vac- 
uum Rabi peaks have the characteristic Lorentzian line 
shape. Their separation and width are given by 2^o and 
(7i+27(^ + /^)/2, respectively. For our specific sample, the 
splitting is observed to exceed 260 linewidths, shown in 
Fig. [2^. Strikingly, when increasing the drive power be- 
yond the linear-response regime, the shape of the trans- 
mission curve changes drastically as shown in Figs. [2]3-f. 
Each vacuum Rabi peak develops a central dip and even- 
tually supersplits into a doublet of peaks. We find excel- 
lent agreement between our prediction based on Eqs. ([2j) 
and ([3]) and the measured transmission. 

The physics of this supersplitting can be elucidated 
within a reduced two-level model. This model takes 
into account only the Jaynes-Cummings ground state 
and either the symmetric or antisymmetric superposition 
of transmon and photon excitation. These two states, 
li) = |0,0) and IT) = (|1,0) ± |0,1))/a/2, form an ef- 
fective two-level system composed of both transmon and 
cavity degrees of freedom. Within the effective two-level 
subspace, the photon operators are mapped^S. to Pauli 
operators a -^ E_/V^, a^ -^ E+/\/2, so that the mi- 
crowave tone acts as a drive on the effective two-level 
system, i.e. 



HA nn 



(4) 



a scenario that Carmichael et alM'^ have referred to 
as "dressing of dressed states." The hamiltonian J^eff 
refers to the frame rotating at the drive frequency, 
A = ciJoi ± ^0 — ^d is the detuning between drive and 
one vacuum Rabi peak, and Q = >/2^ is the effective 
drive strength. With the notable exception of the re- 
cent work by I. Schuster et al.-^, previous investigations 
were primarily concerned with effects on photon corre- 
lations and fluorescence, as observed in photon-counting 
measurement^ii^. According to the operator mapping, 
photon counting can be related to the measurement of 
(S^), whereas detection of the heterodyne transmission 
amplitude A corresponds to |(5]-)|. As a result, hetero- 
dyne detection fundamentally differs from photon count- 
ing and the vacuum Rabi supersplitting is a characteristic 
of heterodyne detection only. 

After restricting the master equation ([2]) to the two- 
level subspace, the system evolution can be expressed in 
terms of simple Bloch equations for the three components 
of the reduced density matrix p = (1 + xT^x + V^y + 

X = -x/T^ -Ay, y = Ax- y/T^ - Qz, 

z = Qy-{z^l)/Ti. (5) 

(An intuitive approach avoiding the Bloch equations 
is discussed in the Supplementary Information, Discus- 
sion.) Here, T{ and Tj are the effective relaxation and 



dephasing times, which are related to 71, 7^,, and k via 
Ti~^ = (71 + k)/2, and T^"^ = (71 + 27^ + «;)/4. Prom 
the steady-state solution of the Bloch equations we find 
the transmission amplitude 



A = 



T^n 



K(A2Tf 



l)/2 



A2r^2 



■ rmn^ + 1 



(6) 



This expression correctly describes the crossover from lin- 
ear response at small driving strength, ft <^ {T[T2)~^^'^ , 
producing a Lorentzian of width 2T2 , to the dou- 
blet structure observed for strong driving. Speciflcally, 
as the drive power is increased, the response saturates 
and the peak broadens, until at 1] = {T[T2)~^^'^ the 
peak undergoes supersplitting with peak-peak separa- 
tion 2T^"V^i^2^ - 1- The fact that we use het- 
erodyne detection is indeed crucial for the supersplit- 
ting. It is easy to verify within the two-level approxima- 
tion that photon-counting always results in a Lorentzian. 
For photon counting, probing beyond the linear-response 
regime merely results in additional power-broadening; 
speciflcally, the width of the Lorentzian is given by 

In Figs. [2]3-e, the analytical expression (|6]) is plotted 
for comparison with the full numerical results and the 
experimental data. We flnd good agreement for low to 
moderate drive power, conflrming that the supersplitting 
can be attributed to driving the vacuum Rabi transition 
into saturation while measuring the transmission with 
the heterodyne technique. For higher drive power a left- 
right asymmetry appears in the true transmission spec- 
trum, which is not reproduced by Eq. (|6]), and which is 
partly due to the influence of levels beyond the two-level 
approximation. 

Higher levels of the extended Jaynes-Cummings hamil- 
tonian become increasingly important when the drive 
power is raised further. Figure [3] shows the emergence 
of additional peaks in the transmission spectrum. Each 
of the peaks can be uniquely identifled with a multi- 
photon transition from the ground state to an excited 
Jaynes-Cummings state. For simplicity, we consider 
the situation where the anharmonicity a and the cou- 
pling strength qq are sufficiently different that mixing 
between higher transmon levels and the regular Jaynes- 
Cummings states |n, ±) = (|n + 1, 0) ± |n, 1))/V^ is min- 
imal for the low-excitation subspaces. Accordingly, the 
experiments are carried out using a different transmon 
qubit in the same sample, with a smaller coupling of 
^o/tt = 94.4 MHz (see Methods). In this case, the n- 
photon transitions to the n-excitation subspace occur at 
frequencies En±/nh = (a;ri^~"'^'^^^o)/27r, and thus reveal 
the nonlinearity of the Jaynes-Cummings ladder. The 
detailed comparison between experimental data and nu- 
merical simulation in Fig. [3]3-e shows superb agreement 
down to the narrowest features observed. 

The possibility of multiphoton transitions at suffi- 
ciently large drive powers also affects the shape of the 
vacuum Rabi splitting when tuning the qubit frequency 
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Figure 3: Emergence of ^/n peaks under strong driving 
of the vacuum Rabi transition, a, The extended Jaynes- 
Cummings energy spectrum. Ah levels are shown to scale 
in the left part of the diagram: black lines represent levels 

|n, zb) ~ (|n + 1, 0) ± |n, 1))/a/2 with only small contributions 
from higher (j > 1) transmon states; grey lines represent lev- 
els with large contributions from higher transmon states. In 
the right part of the diagram, the ^Jn scaling of the split- 
ting between the |n, ±) states is exaggerated for clarity, and 
the transitions observed in plots b— e are indicated at the 
x-coordinate En±/nh of their n-photon transition frequency 
from the ground state, b. Measured transmission (A^, hetero- 
dyne amplitude squared) in colour scale as a function of drive 
frequency and power. The multiphoton transitions shown in 
a are observed at their calculated positions, c— e. Examples of 
cuts for constant power, at the values indicated in b (results 
from the master equation (|2]) in black; experimental results 
in red), demonstrating excellent agreement between theory 
and experiment, which is reinforced in the enlarged insets. 
Good agreement is found over the full range in drive power 
from —45 dB to +3 dB, for a single set of parameters (also see 
Supplementary Movie 1). 
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Figure 4: Qubit— cavity avoided crossing at high drive 
power. Transmission measurement when tuning the qubit 
frequency through resonance for a drive power of +ldB. 
a. Measured transmission as a function of drive frequency and 
magnetic field. As the field is increased, the qubit frequency 
is tuned through resonance with the cavity, and anticrossing 
behaviour is observed. The multiphoton transitions shown in 
Fig. [3^ are visible. The anomaly at B ^ 15.59 is most likely 
due to the crossing of a higher level of the second qubit present 
in the same cavity, b— c. Example cuts at constant magnetic 
field, at the values indicated in a (master equation results, 
calculated using the same parameters as for Fig. [3l are shown 
in black; measured results in red). (Also see Supplementary 
Movie 2). 



c^oi through resonance with the cavity, shown in Fig. [H 
Instead of the simple avoided crossing commonly ob- 
served at low drive powers^, the presence of multiphoton 
transitions leads to a fan-like structure where individ- 
ual branches can again be identified one-to-one with the 
possible transitions in the Jaynes-Cummings ladder. In 
the experimental data of Fig. [4^, processes up to the 5- 
photon transition are clearly visible. Detailed agreement 
with the theory verifies that the more general situation 
of non-zero detuning between qubit and resonator is cor- 
rectly described by our model. 

In summary, we have shown the supersplitting and 
the emergence of multiphoton transitions in the nonlin- 
ear response of the vacuum Rabi resonance, presenting 
both theoretical predictions and their experimental ver- 
ification. This enables the observation of the character- 



istic -Jn nonlinearities in the Jaynes-Cummings ladder 
in a heterodyne transmission measurement, and is di- 
rect evidence for the strong couphng between the su- 
perconducting qubit and the quantized microwave field. 
The high precision of agreement between predictions and 
experiment demonstrates the robustness of the Jaynes- 
Cummings physics in a circuit QED system. This opens 
up pathways for further investigations in quantum optics 
and quantum information, such as generating number- 
squeezed states or employing multi- level quantum logic. 



METHODS 

THEORY 

The modelling of the transmon follows ref. |20|. While 
some expressions given in the main text are asymptotic 
results valid for Ej / Eq ^ 1, our calculations are based 
on a full diagonalization of the Cooper pair box Hamil- 
tonian. 

The description of the transmon-cavity system in 
terms of a master equation requires a model for the re- 
laxation and dephasing of higher transmon levels. As de- 
tailed studies of the microscopic origin of the dominant 
relaxation and dephasing channels are still outstanding, 
we have chosen plausible superoperators for our master 
equation ([2|). Assuming that relaxation of higher trans- 
mon levels may arise due to a coupling of external de- 
grees of freedom to the charge on the superconducting 
island, we take the relative strengths of relaxation to 
be related to the coupling parameters as aj = gj/go- 
Dephasing of higher levels is likely to be due to charge 
noise. Denoting the charge dispersion^^ of level j by 
Cj = (jOjijig = {))—Ui)j{ng = 1/2), we therefore take the rel- 
ative dephasing strengths to be Pj = 2ej/{ei — cq). (The 
normalisation of aj and pj is chosen to allow the usual 
interpretation of 71 and j^ as relaxation and dephasing 
rates for the first two levels of the transmon spectrum.) 
In fact, the pure dephasing rate is sufficiently small for 
our qubits that we set 7^^ = 0. A comparison of addi- 
tional simulations with auxiliary experimental results at 
increased temperatures allows us to place an approximate 
upper bound of 0.003 on the number of thermal photons 
in the cavity. 

For the steady-state solution of Eq. ([2]), the Hilbert 
space is truncated to a subspace with maximum num- 
ber of excitations N, using the projector P/v = 
Zlo<n+j<Ar l^'i) (^'il- I^ o^^ simulations, we keep up 



to A^ = 7 excitations. To reach agreement with the 
experimentally measured signal for the strongest drive 
powers, it is necessary to account for a small amount 
(~ —55 dB) of leakage of the drive past the cavity. In 
addition, there is a small bias introduced by measuring 
the transmission as the square of the / and Q quadra- 
tures, each of which is subject to noise. Accordingly, the 
quantity that corresponds to the experimental signal is 
A^ = \2^tT{psa) + 6^p + 2cr^, where b describes the 
leakage of the drive bypassing the cavity, and a^ is the 
measurement noise in each of the / and Q channels. 

Fits are obtained by minimizing the mean squared 
deviation between experiment and calculation over the 
full power range, with fit parameters being b and the 
two scaling factors describing the signal attenuation and 
amplification for input and output signals. To obtain 
optimal agreement, we also make small adjustments to 
the system parameters 71, /^, go^ cjr, ^01 7 and Ec- These 
parameters are confined by separate experiments to 
narrow ranges, and all values used in fits are consistent 
within the experimental uncertainties. Once obtained, 
the same set of parameters was used in generating Figs. 
[3l m and the Supplementary Movies. 

EXPERIMENT 

Measurements are performed in a dilution refrigerator 
at 15 mK. The sample consists of two superconduct- 
ing transmon qubit^i^S, coupled to an on-chip copla- 
nar waveguide (CPW) cavity. Fabrication of the sam- 
ple followed the description given in ref. 10. The CPW 
resonator has a half-wavelength resonant frequency of 
ujr/27r = 6.92 GHz and a photon decay rate of /T:/27r = 
300 kHz. Transmission measurements are performed us- 
ing a heterodyne detection scheme. The transmitted RF 
voltage signal through the cavity is mixed down to a 
1 MHz carrier signal, and then digitally mixed down to dc 
to obtain the transmitted voltage amplitude as a function 
of frequency. The vacuum Rabi coupling strengths for 
the two qubits are obtained as go/7r = 347 MHz (qubit 1) 
and ^o/tt = 94.4 MHz (qubit 2). Time domain measure- 
ments of the qubits show that they are Purcell-limited 
and completely homogenously broadened at their flux 
sweet spots^^. The coherence times are Ti = 1.7 /is and 
T2 = 0.7 /is (qubit 1, away from the flux sweet spot) and 
Ti = 1.4 /is and T2 = 2.8 /is (qubit 2, at flux sweet spot). 
The charging energies of the two qubits are measured to 
be Ec/h = 400 MHz and Ec/h = 340 MHz. (See Sup- 
plementary Information, Methods, for further details.) 
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